Bragg spectroscopy of a strongly interacting ^^Rb Bose-Einstein condensate 
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We report on measurements of the excitation spectrum of a strongly interacting Bose-Einstein 
condensate (BEC). A magnetic-field Feshbach resonance is used to tune atom-atom interactions in 
the condensate and to reach a regime where quantum depletion and beyond mean-field corrections 
to the condensate chemical potential are significant. We use two-photon Bragg spectroscopy to 
probe the condensate excitation spectrum; our results demonstrate the onset of beyond mean-field 
effects in a gaseous BEC. 



The concept of an interacting but dilute Bose gas 
was originally developed fifty years ago as a theoreti- 
cally tractable surrogate for superfluid liquid helium. For 
some years after the eventual experimental realization of 
dilute-gas BEC, experiments were performed mainly in 
the extreme dilute limit, in which atom-atom correlations 
were of negligible significance. Such correlations again 
assumed a central role, however, in the 2002 experiments 
on a Mott state for bosons in an optical lattice [1] and on 
atom- molecule coherence near a Feshbach resonance 
Atom-atom correlations are also central to the current 
hot-topic field of resonant fermionic condensates [3, 4]. 

In this paper we describe an experimental study of 
elementary excitations in a system which barkens back 
to a previous century, in that, like liquid helium, it is 
a strongly interacting, bulk, bosonic superfluid. Un- 
like liquid helium, our gas of Bose-condensed ^^Rb has 
the modern virtue of Feshbach-tunable interactions well- 
described by a scattering length a that is much larger 
than the reach of the actual interatomic potential. Our 
primary tool for characterizing the sample is Bragg spec- 
troscopy [1,0]. 

The convenience of a Feshbach-tunable scattering 
length in a bosonic system comes at a cost: at high val- 
ues of na^ the gas decays relatively quickly. For this 
reason, we are not able to use, for example, shifts in 
breathing-mode frequencies [3l that have been so fruitful 
in the resonant fermion systems [1, Q . These measure- 
ments require that the gas survive at least as long as 
it takes to traverse the sample at the speed of sound. 
Bragg scattering, on the other hand, is essentially a local 
measurement. The magnetic field can be ramped to tune 
the sample to a large scattering length, and the Bragg 
pulse can be applied on a time scale faster than that 
required for global equilibrium but still consistent with 
local many-body quasi-cquilibrium. 

What should we expect for the energy, fuAj{k), of an ex- 
citation with momentum Hk above a Bose-Einstein con- 
densate of homogeneous density n? In the absence of 
any interactions, the energy of the excitation is simply 
h?k'^ / (2m), where m is the mass of an atom. In the sim- 



plest description of interactions, atoms in the condensate 
feel a chemical potential, /i — 'inh'^na/m jT^]. On the 
other hand, a high momentum excitation has equal di- 
rect and exchange interactions, of magnitude fi each, with 
the low momentum atoms. Thus, the energy required to 
promote an atom out of the condensate, less the bare 
kinetic energy, is 
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The simple model leading to this energy shift, linear in 
n and a (see Fig. 1), relies implicitly on the smallness of 
three dimensionless parameters, which are the following: 



1. VSTrna^ ^ 1, so that interactions are amenable to a 
mean- field treatment. 

2. ka ^ 1, so that the two-body scattering amplitude is 
momentum independent. 

3. l/(fc^) < 1, where ^ = h/{2m^y^'^ is the healing 
length. This limit corresponds to the excitation being 
cleanly in the free-particle regime. For larger l/(fc^), 
the excitation becomes more of a collective object un- 
til, in the opposite limit l/(fc^) 3> 1, the excitation is a 
phonon. The crossover regime was studied experimen- 
tally [III, [lH with the results confirming the original 
theory by Bogoliubov [l3| . 

The primary goal of the work reported in this paper 
is to push into a regime where interactions in a cold 
Bose gas can't be treated as a mean field. For our 
data the next order, Lee Huang Yang (LHY) correction 
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14 [15| is as large as 0.5. 



As we in- 
crease a, the parameters 2 and 3 will also increase. For 
the data presented in this paper l/(fcC) < 0-5 and ka can 
be as large as 0.8. We are not aware of any theoretical 
treatment that simultaneously allows all three parame- 
ters to be nonzero. Understanding the regime in which 
all three parameters are approximately unity is particu- 
larly significant because it corresponds to the intriguing 
roton minimum in the dispersion relation in superfluid 
liquid He [igI , 17 |. The behavior in our gas is in some 
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sense more universal than that in the hquid, because the 
scattering length is much larger than the range of the 
interatomic potential [3 ■ We hope that the results pre- 
sented below, and follow-up work in progress for lower 
values of k, will stimulate interest in the problem. 

In Fig. [2 we plot various calculations good to lowest 
order in combinations of at most two of the three small 
parameters. Beliaev gives a result for ka and l/{k^) be- 
ing nonzero (Eqn. 4.7 in Ref. and numerical in- 
tegration of Eqns. 5.13-5.14 in Ref. [3] (see also Ref. 



[2 d p yields a result for l/{k^) and V Sima^ being nonzero; 
both results are plotted in Fig. 1. We plot also in Fig. 1 
the result of a simple model accounting for finite ka and 
y/Snna^ in the limit of vanishing l/(fc^). Here, Eqn. 1 is 
modified to become 
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inh'^na(k/2) , ^Anfi^na , , 
+ . (2) 



The argument is simply that 1/(A;^) means the exci- 
tation k is very large compared to the momentum content 
of the correlations induced by nonzero v87m(?. Thus, 
for the interaction energy felt by the excitation we ignore 
many-body corrections and simply use S'nh^na{k/2)/m, 
where a{k/2) = -i?e(/(fc/2)) and /(fc/2) is the two-body 
symmetrized scattering amplitude for relative wave vec- 
tor fc/2. For the energy of an atom in the condensate, we 
replace /i with (1 -I- a)Ai:h^na/'m. 

We don't expect Eqn. (2) to be valid in the regime of 
our experiment, because our value of k is insufficiently 
large to ensure the excitations are unaffected by the cor- 
relations induced in the condensate by VoTmc? effects. 
Numerically evaluating the poles of the Green's func- 
tion, Eqns. 5.13-5.14 in Ref. [l3|, shows that our Eqn. 
(2) applies only when l/(fc^) < 0.05, whereas we reach 
l/(fc^) = 0.5. At high momentum (small l/(fc^)), the ef- 
fect of quantum depletion is to shift the dispersion curve 
downwards relative to the condensate energy, while at 
very small momentum, the LHY enhancement to /i in- 
creases the speed of sound and thus increases low-fc val- 
ues of w(fc). The up-shift due to vSTma^ persists even 
to the intermediate momenta of our experiment, which is 
why the dashdotdot curve is above the solid line in Fig. 
1. 
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23. A gas of 85 Rb 
state is first sympa- 



Our experiments are performed using a Rb BEC near 
a Feshbach resonance at 155 G ^ 
atoms in the \F = 2,mp = —2) 
thetically cooled with ^^Rb in a magnetic trap and then 
evaporated directly to ultralow temperature in an optical 
dipole trap [2^. We create a single-species ^^Rh conden- 
sate [2J], with 40,000 atoms and a condensate fraction 
greater than 85%, in a weakly confining optical dipole 
trap at a magnetic field above the Feshbach resonance. 
Curvature of the magnetic field enhances confinement 
along the axial direction of the optical trap. Follow- 
ing evaporative cooling, the optical dipole trap is recom- 
pressed and the final trap has a measured radial (axial) 
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FIG. 1: Theoretical predictions for the peak of the Bragg 
resonance less the bare kinetic energy for an excitation of 
wave vector k = 2 
sity 7.6 X 10^^ cm 
detail in the text and correspond to assuming the following 
quantities are small: ka, VSirna^ and l/(fc^) (green, solid); ka 
and VSTvna^ (blue, dots); ka (magenta, dashdotdot); VSnna''^ 
(red, dashes); and l/(fcC) (black, dashdot). For reference the 
top axis shows the size of the LHY correction, a, relative to 
the condensate chemical potential. 



trap frequency of 27r x 134 Hz (27r x 2.9 Hz), yielding a 
condensate mean density of 2.1 x 10^'^ cm"'^. 

Bragg spectroscopy via stimulated two-photon transi- 
tions provides a direct probe of the condensate excitation 
spectrum. Two counter-propagating, near-resonant laser 
beams are aligned along the long axis of the condensate. 
The momentum imparted to a condensate excitation is 

is the wave 



given by hk = 2hkL where fc^ = =57, 
vector of a beam. The excitation energy is scanned by 
adjusting the frequency difference of the two laser beams. 
The average of the two frequencies is red detuned from 
atomic resonance by 4.2 GHz. The intensity and pulse 
duration of the Bragg beams are chosen so that the frac- 
tion of the condensate excited is less than 10%. 

Just before performing the Bragg spectroscopy, we 
transiently enhance the condensate density by means of 
large amplitude radial and axial breathing modes, which 
we excite by modulating the magnetic field and thus the 
Feshbach-modified scattering length. The rates of the 
ramps are limited so the d/a never exceeds 0.06fi./ (ma^). 
The scattering length is derived from measurements of 
the magnetic field and a previous measurement of the 
85Rb Feshbach resonance [2^. Sychronized with the 
inner turning point of the radial oscillation, we ramp 
the scattering length to the value for a given measure- 
ment and then pulse on the Bragg beams. During the 
pulse, the cloud's inward motion is checked and it be- 
gins to breathe outward. We model the resulting time- 
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dependent condensate density using a variational solu- 
tion to the Gross-Pitaevskii equation [2^ , which predicts 
that the density of the cloud does not change by more 
than 30% during the Bragg pulse. We can meet this 
goal only by using progressively shorter Bragg pulses for 
higher values of desired a. The time- and space-averaged 
density during the pulse is approximately 7.6x 10^'^ cm^'^, 
but this depends weakly on the final value of a. 

After the Bragg pulse, we ramp a to a large, fixed 
value to ensure that the momentum of the excitations 
is spread via collisions 27, [i^ to the entire condensate 
sample. We then infer the total momentum, and thus 
excitation fraction, from the amplitude of the resulting 
axial slosh, measured via an absorption image taken of 
the cloud at a time near its axial turning point. 
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FIG. 2: Typical Bragg spectra at a scattering length of 100 
ao (blue triangles), 585 ao (red circles), and 890 ao (black 
squares) . The excitation fraction is determined from the mea- 
sured momentum transferred to the BEG and plotted as a 
function of the frequency difference between the two Bragg 
beams. Lines are fits of the data as described in the text. 
Mean-field theory predicts a continuous increase in the line 
shift with increasing a, however by 890 ao our data display a 
decreasing shift with stronger interactions. 



Figure[2]shows measured Bragg spectra for three values 
of a. We fit each Bragg spectrum to an antisymmetric 
function assuming a Gaussian peak and extract a center 
frequency and an RMS width. The Bragg line shift is 
the difference between the fitted center and the ideal gas 
^'^^ - 15.423 kHz. In Fig. EJi we plot our 
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measured line shifts as a function of the scattering length 
a. For a < 300ao (where the predicted LHY correction is 
already a 10% effect), the measured line shift (• in Fig. 
[3^) agrees with the simple mean- field result (Eqn. 1). 
However, as the scattering length is increased further, 
the resonance line shift deviates significantly from the 
mean-field prediction. The measured line shift reaches 
a maximum near a = 500ao and then decreases as the 
scattering length is increased further. 

At large a we find that our measured line shift exhibits 
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FIG. 3; (a) Bragg line shift and (b) width as a function of 
scattering length. The hollow circles are our observations. 
The solid circles are data corrected for a fitting systematic 
associated with the broad thermal atom background, and the 
error bars represent fit uncertainties. The theory lines and 
LHY correction in (a) are as in Fig. 1 except they are calcu- 
lated for the trapped gas using a local density approximation 
for each of the corresponding data points; the mean BEG den- 



sity ranges from 6.3 x 10" cm ^ to 7.6 x 10 



In (b) 



the solid black circles are the rms width of a gaussian fit to 
the Bragg spectra. Black triangles are from a fit to a convo- 
lution of various contributions to the width calculated under 
the conditions of our measurements. The remaining symbols 
characterize constituent contributions to the convolution in- 
cluding the Lorentzian FWHM width due to collisions (blue 
squares), and the RMS widths of contributions due to the in- 
homogeneous density (red diamonds) and the pulse duration 
(green circles). The largest contribution to the width comes 
from the pulse duration; because the jump to large a initiates 
rapid expansion of the BEG, ever shorter pulses are used to 
obtain the spectra at larger a. 



a systematic dependence on the temperature of the sam- 
ple [2§]. Non-condensed ^^Rb atoms also respond to the 
Bragg pulse, and this causes an observable effect in the 
measured line shift when the spectral width of the con- 
densate response becomes comparable to that of the non- 
condensed atoms (for a > 500ao). We vary the temper- 
ature of the gas to characterize this effect and we apply 
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a small correction to our data to represent the expected 
line shift at zero temperature (• in Fig. [3^). 

The theory curves of Fig. 1 are adapted for direct com- 
parison with our data (Fig 3a) by calculating, at each 
tick mark on the curves, a spatial and temporal average 
over a local density approximation for the correspond- 
ing data point, assuming the density of the condensate is 
given by a time-varying Thomas-Fermi profile. In these 
predictions, we account for observed number losses (typ- 
ically < 30%) during our experiments. The observed loss 
is consistent with the previously observed a-dependent 
three-body recombination rate of ®^Rb (2^ . 

Figure[3|3 shows the measured width of the Bragg peak 
as a function of a. Several effects contribute to the to- 
tal width of the Bragg resonance including the finite du- 
ration of the pulse, the inhomogeneous density of the 
trapped condensate |5|1, collisions between the excitations 
and the condensate [27|, and Doppler broadening [5;]. In 
our case, Doppler broadening is negligible since the ax- 
ial size of the condensate is relatively large. To under- 
stand the total width we convolve the various calculated 
lineshapes of the remaining three effects. We fit a Gaus- 
sian to the convolution (to match the Gaussian fit to our 
data) and the RMS width from this fit is shown in Fig. 
|3Jd (black A). Fig. [Sb also shows the expected contribu- 
tions to the width from each of the three effects. For the 
lineshape due to collisions we expect a Lorentzian with a 
full width at half maximum Si/ = tt- ^^^^ , where a is the 
elastic cross section for collisions between the excitations 
and low momentum atoms. In calculating a we include 
the suppression in the phonon regime predicted by Be- 
liaev [H, The measured Bragg width exceeds the 

predicted width in the strongly interacting regime. How- 
ever, many of the theoretical difficulties in describing the 
line shift apply also to predicting the width that arises 
from inhomogeneous density and excitation lifetime. 

In comparing our measured line shifts with the vari- 
ous theory models (Fig. 3), one should remember that 
none of the models takes into account the fact that the 
inter- atomic potential supports bound states, whose ex- 
istence means that a Bose-condensed gas can be only a 
quasi-equilibrium state for any atomic species. This is 
especially relevant for increasing values of a. A key fu- 
ture goal of our work is to experimentally explore differ- 
ent time-scales for the establishment of local many-body 
quasi-equilibrium, and for longer-term evolution. The ex- 
treme aspect ratio of our sample hastens the loss of den- 
sity that occurs during the expansion caused by the ramp 
to high a. At present, we are required to use short du- 
ration Bragg pulses, which limits our spectral resolution, 
and we are prevented from tracking the time evolution of 
line shifts. An ongoing redesign to a more spherical ge- 
ometry will help. In addition, the Bragg beams are being 
reconfigured to allow access to the low-fc, pure-phonon 
regime for which l/(fc^) ^ 1. 
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